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ABSTRACT 

Since  the  crystal  field  can  be  considered  as  a  perturbation 

of  the  free  ion  spectra  of  actinide  ions  in  certain  crystals,  we  have 

3 

calculated  energy  levels  and  wave  functions  of  the  5f  configuration, 

3+ 

which  is  the  ground  configuration  of  U  ,  in  intermediate  coupling  for 
values  of  y  =  £  /F^  of  7  to  11.  Energy  levels  are  given  in  units  of  F^, 
and  wave  functions  are  expressed  as  linear  combinations  of  Rus  sell- 

5  t  Saunders  functions.  Hydrogenic  and  non-hydrogenic  approximations 

■< 

j  are  compared. 

* 

Reduced  matrix  elements  in  intermediate  coupling  are  given 

so  that  crystal  field  splittings  may  be  calculated.  The  splittings  of 
,  3+ 

the  two  lowest  levels  of  U  in  a  weak  field  with  and  without  J -mixing 
are  compared. 


L 


INTRODUCTION 


In  the  triply  ionized  actinides  the  optically  active  electrons 
n  1  Z 

belong  to  the  5f  configuration  ’  ,  so  comparison  of  their  spectra 

with  those  of  the  triply  ionized  rare  earths,  whose  optically  active 

electrons  belong  to  the  4fn  configurations,  is  of  considerable  interest. 

The  effect  of  the  crystal  field  on  rare  earth  ions  is  small  enough  to 

be  treated  as  a  perturbation  of  the  free  ion  spectra;  this  is  also 

3-9 

true  for  the  actinide  ions,  at  least  in  some  crystals  ,  In  a  free 

ion,  electric  dipole  transitions  within  a  configuration  are  forbidden, 

but  the  odd  parity  part  of  the  crystal  potential  will  mix  configurations 

of  opposite  parity  so  that  "forced"  electric  dipole  transitions  can 

take  place  as  well  as  the  allowed  magnetic  dipole  transitions.  The 

strength  of  the  forced  electric  dipole. transition  depends  inversely 

10  11 

upon  the  energy  differences  between  mixing  configurations  ’  .  In 

crystals  containing  triply  ionized  actinides  total  absorption  begins  at 

about  4000A,  and  it  is  assumed  that  this  is  due  to  the  allowed  5f-6d 

transitions.  (Knowledge  of  the  free  ion  spectra  would  be  useful  here). 

Since  the  analogous  4f-5d  transitions  in  rare  earth  crystals  take 

1  2 

place  at  about  2000  A  ,  more  intense  spectra  can  be  expected  from 

this  effect  in  actinide  crystals  than  in  rare  earth  crystals. 

3+  3 

The  ground  configuration  of  U  is  5f  .  Since  nothing  is 

3+ 

known  experimentally  about  the  free  ion  spectrum  of  U  ,  we  have 

3 

made  calculations  on  the  5f  configuration  to  aid  in  the  interpretation 
of  crystal  spectra. 
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II.  FREE  ION 

The  Hamiltonian  for  the  free  ion  is 

H,  =  H  +  H  +  H 
foes 


where  H  is  the  interaction  of  the  5f  electrons  with  the  core  (constant 
o 

for  the  configuration) 

H  is  the  electrostatic  interaction  between  the  three  5f 
e 

electrons 


Hg  is  the  spin  orbit  interaction. 

Higher  order  effects,  including  spin-spin  and  spin-other -orbit 

interactions  will  be  neglected. 

The  electrostatic  interaction  is 
2 


H 


I  h  i>j  =  *’2’ 


i>j 


where  electrons  in  closed  shells  are  assumed  to  affect  the  energies 

of  the.  5f  electrons  only  by  a  contribution  to  the  central  potential. 

1  13 

We  can  expand  -  in  terms  of  the  Legendre  polynomials 

rij 

cd  n 


r . . 
iJ 


)  —  P  (cos  (1)  ) 

L  n  n  ' 


k=0 


where  r  is  the  smaller  and  r,  the  larger  of  the  distances  of  r.  and 
a  b  °  i 

r j  from  the  nucleus,  and  ou  is  the  angle  between  r^  and  r^.  Then  we 
can  separate  the  radial  and  angular  components  and  write 


i  k  j;  o  5 
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H  =  )  fnF  n  =  2,  4,  6  for  f  electrons 

e  L  n 

n 

Matrix  elements  of  frl  are  essentially  integrals  of  spherical  harmonics 

14 

over  the  angular  wave  function,  and  have  been  evaluated  by  Carlson 

3  ra 

for  f  using  Racah's  method.  The  F  are  integrals  of ^  over  the 

5f  radial  wave  functions  (Slater  integrals).  If  5f  hydrogenic  wave 

15 

functions  are  assumed,  that  is,  in  a  Coulomb  field,  Judd  has 
calculated 

F4/F2  =  0.1422  and  F6/F2  =  0.  0161 

The  spin  orbit  interaction  is 
3 


i=  1 


£,,  the  integral  of  £(r^)  over  the  radial  wave  function,  is  determined 
experimentally  due  to  lack  of  knowledge  of  the  wave  function.  The 
matrix  elements  of  the  angular  part  of  Hg  have  been  calculated  for 
f^  by  Judd  and  Loudon*^.  In  U^+  the  spin  orbit  interaction  is  rather 
large,  therefore  we  had  to  perform  intermediate  coupling  calculations. 
Since  J  is  a  good  quantum  number,  this  involved  diagonalizing 
matrices  for  each  J  value,  l/2  to  17/2.  The  electrostatic  matrices, 
assuming  5f  hydrogenic  wave  functions,  and  spin  orbit  matrices 
have  been  diagonalized  simultaneously,  by  using  a  program  developed 


1 


lrt> 

o 

K 

w 
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by  Wybourne*,  for  various  values  of  y  =  4/F2  exPres®ing  energies 
in  terms  of  and  wave  functions  as  linear  combinations  of  Russell- 
Saunders  states.  Figure  1  shows  the  variation  in  energy  levels  of 

3 

f  with  for  y  =  7  to  11.  The  levels  are  named  by  the  Russell  - 

Saunders  states  into  which  they  go  at  4  =  0.  For  F^  =  200,  y  =  8, 

3+17 

(the  order  of  magnitude  appropriate  for  U  )  ,  it  is  noted  that  the 

-1  3+  3 

configuration  extends  35,  000  cm  .  In  comparison,  the  Nd  4f 

-1 

configuration  extends  67,  000  cm  .  The  41  levels  are  each 
(2J  +  l)-fold  degenerate. 

3+ 

The  use  of  the  5f  hydrogenic  wave  functions  for  U  is 

g 

questionable.  Lammermann  and  Conway  feel  that  the  Coulomb 

3+  5 

approximation  is  satisfactory  for  the  interpretation  of  the  Pu  (5f  ) 
spectrum.  Cohen  has  done  a  relativistic  self-consistent  calculation 
for  the  normal  uranium  ion,  and  Winocur  has  calculated  the 
following  set  of  Slater  integrals  from  Cohen's  eigenfunctions*® 

F4/f2^°.159  F6/F2  =  0.  0204. 

The  discrepancy  between  a  similar  calculation  for  PrIV  and  the 

experimental  values  indicates  that  these  should  be  treated  with 
17 

reserve 

9 

McLaughlin  determined  F.,,  F^,  and  F^  experimentally 

2 

by  a  least  squares  fit  to  the  spectrum  of  UCl^(5f  )  and  found  the 
ratios  to  be 

*  Private  communication. 


! 


-5- 


F  ./F,  =  0.  1468 


F6/F2  =  0.  0219 


We  have  calculated  energy  levels  and  wave  functions  of  f  using 

these  last  ratios  with  y  =  8.0,  and  the  energy  levels  are  compared 

with  hydrogenic  ones  in  Figure  2.  F2  has  been  chosen  to  put  the 
4  -1 

*ll/2  leve*  at  4500  cm  from  the  ground  state  where  it  has  been 


19  20  4 

observed  ’  It  is  apparent  that  the  position  of  the  F 


3/2 


level 


relative  to  the  1^3/2  *8  depressed  in  the  non -hydrogenic  case, 
and  above  10,  000  cm  ,  where  the  energy  levels  become  more 
dense,  the  grouping  is  much  different  in  the  two  cases. 
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III.  CRYSTAL  FIELD 

We  are  interested  in  crystals  in  which  the  splitting  of  the  J 

manifold  is  small  compared  to  the  spin-orbit  splitting,  so  in  a  first 

approximation  the  crystal  field  interaction  can  be  treated  as  a 

perturbation  of  the  free  ion  energy  levels.  However,  since  the  5f 

22 

wave  functions  are  more  extensive  than  the  4f  wave  functions  , 

somewhat  larger  crystal  field  effects  can  be  expected  in  the  actinides 

than  in  the  rare  earths.  Each  J  level  of  an  ion  with  an  odd  number 

of  electrons  will  split  in  a  crystal  field  into  ( 2 J  +  l)/2  components 

at  the  most  due  to  Kramers  degeneracy. 

The  Hamiltonian  for  the  crystal  field  interaction,  assuming 
23 

a  purely  electrostatic  field  is 

H  =  Z  (r^)  V^(8,  4>)  |q|  <  k,  k  =  2,  4,  6  for  f  electrons. 


The  number  of  terms  in  the  expansion  depends  on  the  symmetry  of 

the  ion  site  in  the  crystal.  The  coefficients  A^  depend  on  the  lattice, 

and,  although  it  is  theoretically  possible  to  calculate  them  if  the 

positions  and  charge  distributions  of  the  ions  in  the  crystal  are 

known,  results  have  not  been  satisfactory.  Therefore,  A^  (  r^) 

are  customarily  treated  as  experimentally  determined  parameters. 

The  are  normalized  associate  spherical  harmonics  which  transform 

in  the  same  manner  as  the  unit  tensor  operators  and  are  related 

to  them  by  constant  factors.  We  may  use  either  set  for  calculations. 

Matrix  elements  of  the  crystal  potential,  for  LS  coupling  then,  can 

15 

be  obtained  by  using  the  relation 


I 
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,—<* 


<aLSJM|u^|a'L'S'J'M')=  6{SS'H-l)J'M(.MqM,)(aLSJ  |  |uk|  |a'L'S'J'>, 


J  kJ' 


where  the  reduced  matrix  element  is 


<pLSJ  ||uK||a'L’S'J')=(-l) 
J  kJ 


S+k+L'+J[(2J+1)(2J'  +  l)]  /  {^L,^}(aLS||Uk||a'L'S). 

L  J  S. 


o  V-/ 

3-j  symbols  (  an<^  6-j  symbols  '  jt  Li  ^  j  are  vector  cotipling 


coefficients  and  are  available  in  tabular  form  Judd  has  tabulated 
the  double  reduced  matrix  elements  (  aLS  ||  U^Ha'L'S)  for  f^.  If  J 
mixing  is  neglected,  the  matrix  elements  may  be  more  simply 
calculated  from 


(aLSJM |  | a' L'SJM' )  =  <  a  LSJ  1 1  Vk  |  |a 'L'SJ  >  fks£(M), 

where  S^(M)  are  the  Stevens  coefficients  and  the  fk>  which  depend 
only  on  J,  are  proportionality  constants  chosen  to  make  the  S° 

K 

integers . 

3+ 

U  states  can  be  represented  by  linear  combinations  of 

Russ  ell -Saunders  states  with  coefficients  depending  on  y  ;  therefore 

the  crystal  field  matrix  elements  will  also  depend  on  y  ,  and  we  have 

obtained  them  for  A  J  =  0  by  the  same  transformation  which  diagonalized 
3 

the  5f  energy  matrix. 

Table  I  presents  the  results  of  calculations  we  have  made 
3 

on  the  5f  configuration  using  hydrogenic  wave  functions.  For  each 
J  value  for  y=  7.0  to  11.0,  energies  in  units  of  F^  are  given  in  the 
first  row.  In  the  columns  below  each  energy  level  a  re  the  coefficients 
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of  the  LS  wave  functions  for  that  level.  Reduced  matrix  elements 

in  intermediate  coupling  are  given  in  terms  of  (  aLSJ  1 1  |  |aL'SJ  ) 

3 

f^  X  10  in  order  that  they  may  be  used  with  the  Stevens  coefficients 
26  27 

given  in  Table  II  ’  The  last  row  of  each  column  gives  the 
Lande*  g -value  of  the  level  in  intermediate  coupling. 

As  an  example  of  the  use  of  the  tables,  the  splitting  of 
a  level  in  symmetry  for  y  =  8  is  ( -0 . 00937)  S^( M)  + 

(-0.  02391)  S°(M)  or 

(-0.  00937)X5XA°+(-0.  02  391)XlXA°=-0.  04685 A^-  0. 02391A°,  m=5/2 
-(0.00937)X(-1)XA°-(0.02  391)X(-3)XA°=0.  00937A°+0.  07173A®,  m=3/2 
-(0.00937)X(-4)XA°-(0.02391)X2  A°=0.  03748  A°  -0 . 04782  A°,  m=l/z  . 


In  a  crystal  field  J  is  only  approximately  a  good  quantum 

number,  so  that  J  levels  will  also  be  mixed.  A  complete  J  mixing 

calculation  would  involve  diagonalizing  three  60-  and  61 -dimensional 

matrices,  and  does  not  seem  to  be  worthwhile  until  more  is  known 

about  the  wave  functions.  However,  in  order  to  get  an  idea  of  the 

4  4 

magnitude  involved,  J  mixing  of  the  Ig/2  ant^  *ll/2  s*-a^e9  can  be 
calculated  by  second  order  perturbation  theory.  (For  higher  levels 

which  become  more  nearly  degenerate  this  cannot  be  done). 

4  4 

We  have  calculated  the  splittings  of  the  *9/2  and 
states  using  5f  hydrogenic  wave  functions  for  y  =  8,  0  and  the  fictitious 
crystal  field  parameters 
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A°<  r2)  =  170 


2 

Aj  <  r4)  -  -80 


A^(rS=  -100 

A^  (r6>  =  1100 


4  3+ 

which  are  close  to  those  found  by  Gruber  for  Am  in  LaCl^. 

4  4  -1 

The  effect  of  the  I state  on  the  1^ state  4500  cm  below 

it  is  to  depress  the  ground  Stark  level  by  27  cm  and  to  increase 


slightly  the  overall  splitting.  The  calculated  splittings  of  the  I 
4 

and  Ig/2  levels  using  the  above  parameters,  with  and  without 
J  mixing,  are  shown  in  Figure  3. 


ll/2 


JL  IV 
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CONCLUSION 


Experimentally,  a  complication  occurs  in  uranium  doped 

4+ 

crystals  in  that  uranium  readily  ionizes  to  U  and  less  easily  to 

,r3+ 

ing  U  s 
3+  28, 29 


U^+  so  that  crystals  nominally  containing  U^+  sometimes  contain 


both  ions,  as  was  the  case  for  CaF^iU 


The  ground 


4+  2 

configuration  of  U  is  5f  .  In  Figure  4  we  compare  McLaughlin's 

4+  3 

energy  level  diagram  for  U  with  our  calculated  one  for  5f  for 

reasonable  values  of  F^  and  y,  and  find  that  the  energy  levels 

are  often  very  close,  especially  when  we  take  into  account  Stark 

splitting,  which  could  spread  each  level  by  several  hundred  cm 

However,  a  1.88p  absorption  line,  described  as  characteristic 
3+  28  29  3 

of  U  ’  does  not  fit  into  the  5f  energy  level  scheme.  Either 

hydrogenic  wave  functions  are  decidedly  in  error  or  this  1.88p  line 

4+  3 

belongs  to  U  ,  where  it  could  fit  into  the  state. 

It  is  difficult  to  see  how  these  levels  can  be  definitely 

assigned  until  the  free  ion  spectra  have  been  analyzed.  Also  it 

is  quite  desirable  to  know  the  position  of  the  6d  configuration  in 

order  to  estimate  line  strengths.  Therefore,  a  program  has  been 

initiated  to  analyze  the  free  ion  spectra  of  ionized  uranium.  We 

3+ 

hope  to  start  with  UIV  (U  )  and  to  go  on  to  UIII  and  UV.  This 
work  will  be  done  in  cooperation  with  the  Physics  Department  of 
The  Johns  Hopkins  University. 


f 

i 


} 
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FOR  y=  8. 

(a)  F?  =  211  F,  -  30.  004  F,  =  3.  3971  (hydrogenic) 

(b)  F  =  213  F?  =  31 . 098  .  F,  =  4.  6647  (McLaughlin's  ratios 

L  4  from  U  ), 


I 


-80  A,  (r  >=1100 
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6.904  105 
5.765  875 
26.75?  484 
117.166  548 
58.584  024 

1.857  559 
16.911  555 
84.095  184 
29.752  158 

7.946  248 
51.497  575 
20.752  510 

85.551  15* 
12.401  997 

5.546  542 

(b)  S66 


TABLE  II.  STEVENS  COEFFICIENTS 


ADDITIONAL  TABLES  OF  S^  MAY  BE  FOUND  IN  REFERENCE  2. 
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